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Abstract 



The entropy is one of the most applicable uncertainty measures in many statistical and en- 
gineering problems. In statistical literature, the entropy is used in calculation of the Kullback- 
Leibler (KL) information which is a powerful mean for performing goodness of fit tests. Ranked 
Set Sampling (RSS) seems to provide improved estimators of many parameters of the popu- 
lation in the huge studied problems in the literature. It is developed for situations where the 
variable of interest is difficult or expensive to measure, but where ranking in small sub-samples 
is easy. In This paper, we introduced two estimators for the entropy and compare them with 
each other and the estimator of the entropy in Simple Random Sampling (SRS) in the sense 
of bias and Root of Mean Square Errors (RMSE). It is observed that the RSS scheme would 
improve this estimator. The best estimator of the entropy is used along with the estimator of 
the mean and two biased and unbiased estimators of variance based on RSS scheme, to esti- 
mate the KL information and perform goodness of fit tests for exponent iality and normality. 
The desired critical values and powers are calculated. It is also observed that RSS estimators 
would increase powers. 

Keywords: Ordered Ranked set sampling; Judgement ranking; Order statistic; Information theory; 
Exponential; Normal; Uniform 

1 Introduction 

Suppose a continuous random variable X has cumulative distribution function (cdf) F(x) and a 
probability density function (pdf) f(x). The differential entropy H(f) of the random variable X 
is defined to be 
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The entropy is one of the most applicable uncertainty measures in many statistical and engineering 
problems. In statistical literature, the entropy is used in calculation of the Kullback-Leibler (KL) 
information which is a powerful mean for performing goodness of fit tests. The Kullback-Leibler 
(K-L) information of f(x) against /o(a;) is defined in [7] to be 



Since /(/; /o) has the property that I(f;fo) > 0, and the equality holds if only if / = /o, the 
estimate of the K-L information has also been considered as a goodness of fit test statistic by 
some authors including [2] and [5]. It has been shown in the aforementioned papers that the 
test statistics based on the K-L information perform very well for testing exponentiality [5] as 
compared, in terms of power, with some leading test statistics. 

Ranked Set Sampling (RSS) has been developed by Mclntyre (1952). This method is applied 
for situations in which measuring a variable is costly or difficult, but where ranking in small subsets 
is easy. In this method, we first subdivide a sample of size n — k 2 randomly into k subsamples 
of size k, rank each subsample visually or using any simple or cheap method and then in the r th 
subsample, measure and record only the unit of rank r which is denoted by X^l (r — 1, . . . , k). 

(r) a 

Since the subsamples are independent, X^.^'s are independent random variables. Also the marginal 
distribution of X^J. is the same as that of r* n order statistic from a sample of size k of X, i.e. 
X r:k . As it was proved by Mclntyre, mean of this sample is an unbiased estimator of the mean of 
Y with an efficiency slightly less than \{k + 1), relative to the mean of a Simple Random Sample 
(SRS) of size k. Thus "ranked set sampling should be useful when the quantification of an clement 
is difficult but the elements of a set are easily drawn and ranked by judgment." (Dell and Clutter 



This method was also extended to estimating variance (Stokes 1980a), correlation coefficient 
(Stokes 1980b) and the situations in which the sample is subdivided into subsamples of different 
sizes. 

In This paper, we introduced two estimators for the entropy and compare them with each other 
and the estimator of the entropy in Simple Random Sampling (SRS) in the sense of bias and Root 
of Mean Square Errors (RMSE). It is observed that the RSS scheme would improve this estimator. 
The best estimator of the entropy is used along with the estimator of the mean and two biased and 
unbiased estimators of variance based on RSS scheme, to estimate the KL information and perform 
goodness of fit tests for exponentiality and normality. The desired critical values and powers are 
calculated. It is also observed that RSS estimators would increase powers. 

2 Entropy estimation 

The nonparametric estimation of the entropy 
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1972). 
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Table 1: Simulated Minimum RMSE (MRMSE) and Minimum Absolute Bias (MAB) of H^ n and 
R mn and optimal to for k = 10 and three distributions with different values of r. 







2 


r 


3 








ttI 

mn 


mn 


ttI 

mn 


tt2 

mn 


U(0,1) 


MRMSE (optimal m*) 


0.062(8) 


0.081(5) 


0.045(11-13) 


0.073(5) 




MAB (optimal m*) 


0.030(10) 


0.047(5) 


0.021(15) 


0.048(5) 


e(l) 


MRMSE (optimal to*) 


0.157(5) 


0.168(4) 


0.125(6) 


0.140(4) 




MAB (optimal to*) 


0.001(6) 


0.0137(5) 


0.004(7) 


0.014(5) 


N(0,1) 


MRMSE (optimal to*) 


0.184(5,10) 


0.246(5) 


0.138(7,8) 


0.233(5) 




MAB (optimal to*) 


0.113(10) 


0.205(5) 


0.062(12) 


0.206(5) 



m = l(l)ft/2 for H? nn and m = l(l)rfc/2 for Hi 



An estimate of ([3]) can be constructed by replacing the distribution function F by the empirical 
distribution F n . The derivative of F~ l (i/n) is estimated by (xi+ w:n — Xi- w:n )n/(2w). The estimate 
of H is then 

1 n n 

H(m,n) = - ^ log {^^{x l+m , n - x^ m:n )j , (4) 

i— 1 

where the window size m is a positive integer, which is less than n/2, and Xi :n = x\ :n for i < 1, 
and a^ : „ = a;„ : „ for i > n. 

Ebrahimi et al. (1994) proposed a modified sample entropy as 

n 

H c (n,m) = n^^log (X (l+m) - X (l _ m) ) (5) 

CiTfl 

where 

1 + i=i if 1< i < m 

m — — 

c.j = ^ 2 if to + 1 < i < n — to 

1 + 2=i tf n -m + l<i<n 

m — — 

To estimate the entropy in RSS scheme, we may note that the estimator of F _1 (i/n) must be 
positive for log function to be well-defined. So we have to order the ranked set sample. There are 
two ways to order this sample. First way is to order each replication, derive the estimator and 
then take the average as the main estimator. The second way is to order the whole sample of size 
rk. This two methods yield two estimators as follows 

1 n n 

H ^n = - l0 § TZZ ( X li+m] ~ X[i-m] ) (6) 
i=l 1 

and 

1 r k k 

3=1 i = l 1 
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Table 2: Monte Carlo biases and RMSE for H^n in three distributions for n = 10, 20 



U(0,1) 



e(l) 



N(0,1) 



n 


m 


SRS 




RSS 




SRS 




RSS 




SRS 




RSS 




Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


10 


1 


-0.381 


0.451 


-0.259 


0.326 


-0.392 


0.561 


-0.298 


0.398 


-0.452 


0.458 


-0.342 


0.428 




2 


-0.222 


0.293 


-0.108 


0.168 


-0.222 


0.436 


-0.142 


0.266 


-0.342 


0.441 


-0.227 


0.311 




3 


-0.159 


0.228 


-0.070 


0.124 


-0.174 


0.405 


-0.078 


0.241 


-0.301 


0.408 


-0.207 


0.289 




4 


-0.140 


0.224 


-0.056 


0.110 


-0.114 


0.382 


-0.031 


0.236 


-0.305 


0.394 


-0.209 


0.285 




5 


-0.131 


0.212 


-0.050 


0.107 


-0.064 


0.371 


0.012 


0.249 


-0.289 


0.389 


-0.204 


0.279 


20* 


1 


-0.328 


0.358 


-0.274 


0.302 


-0.335 


0.424 


-0.290 


0.340 


-0.373 


0.427 


-0.313 


0.358 




2 


-0.176 


0.203 


-0.121 


0.147 


-0.179 


0.299 


-0.139 


0.206 


-0.221 


0.288 


-0.178 


0.232 




3 


-0.125 


0.155 


-0.076 


0.103 


-0.151 


0.280 


-0.083 


0.169 


-0.179 


0.252 


-0.141 


0.200 




4 


-0.104 


0.134 


-0.056 


0.084 


-0.098 


0.264 


-0.052 


0.161 


-0.176 


0.255 


-0.124 


0.189 




5 


-0.088 


0.119 


-0.046 


0.074 


-0.062 


0.253 


-0.024 


0.157 


-0.167 


0.245 


-0.117 


0.184 




6 


-0.079 


0.117 


-0.040 


0.067 


-0.047 


0.244 


0.001 


0.165 


-0.156 


0.232 


-0.116 


0.185 




7 


-0.076 


0.111 


-0.035 


0.063 


-0.020 


0.263 


0.025 


0.173 


-0.150 


0.231 


-0.116 


0.185 




8 


-0.068 


0.109 


-0.034 


0.062 


0.010 


0.264 


0.051 


0.188 


-0.157 


0.239 


-0.114 


0.185 




9 


-0.064 


0.108 


-0.032 


0.063 


0.032 


0.260 


0.078 


0.203 


-0.158 


0.241 


-0.116 


0.188 




10 


-0.061 


0.106 


-0.030 


0.063 


0.044 


0.268 


0.102 


0.225 


-0.152 


0.234 


-0.113 


0.184 



¥ n = lOr cases are observed by RSS scheme with 10 samples and r replication. 



Tabic 3: Monte Carlo biases and RMSE for H mn in three distributions for n — 30 



U(0,1) 



e(l) 



N(0,1) 



n to 


SRS 




RSS 




SRS 




RSS 




SRS 




RSS 




Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


Bias 


RMSE 


30* 1 


-0.312 


0.337 


-0.273 


0.290 


-0.293 


0.363 


-0.286 


0.318 


-0.328 


0.370 


-0.300 


0.331 


2 


-0.158 


0.174 


-0.125 


0.141 


-0.156 


0.248 


-0.136 


0.181 


-0.198 


0.245 


-0.160 


0.199 


3 


-0.110 


0.129 


-0.080 


0.096 


-0.115 


0.219 


-0.084 


0.145 


-0.158 


0.209 


-0.118 


0.164 


4 


-0.090 


0.110 


-0.058 


0.075 


-0.078 


0.206 


-0.052 


0.129 


-0.135 


0.196 


-0.099 


0.150 


5 


-0.071 


0.092 


-0.046 


0.065 


-0.059 


0.202 


-0.034 


0.126 


-0.113 


0.184 


-0.088 


0.144 


6 


-0.069 


0.090 


-0.039 


0.058 


-0.045 


0.198 


-0.013 


0.125 


-0.106 


0.181 


-0.082 


0.142 


7 


-0.061 


0.083 


-0.034 


0.054 


-0.023 


0.206 


0.004 


0.129 


-0.098 


0.174 


-0.073 


0.138 


8 


-0.056 


0.079 


-0.030 


0.050 


-0.007 


0.194 


0.021 


0.135 


-0.106 


0.175 


-0.069 


0.138 


9 


-0.052 


0.076 


-0.028 


0.048 


0.015 


0.192 


0.039 


0.145 


-0.086 


0.174 


-0.067 


0.140 


10 


-0.050 


0.075 


-0.027 


0.046 


0.027 


0.195 


0.057 


0.155 


-0.091 


0.175 


-0.067 


0.141 


11 


-0.048 


0.075 


-0.025 


0.045 


0.050 


0.211 


0.073 


0.167 


-0.091 


0.171 


-0.065 


0.141 


12 


-0.041 


0.071 


-0.024 


0.045 


0.075 


0.225 


0.098 


0.185 


-0.090 


0.171 


-0.062 


0.141 


13 


-0.042 


0.074 


-0.023 


0.045 


0.089 


0.231 


0.117 


0.201 


-0.089 


0.175 


-0.065 


0.144 


14 


-0.043 


0.073 


-0.022 


0.046 


0.100 


0.248 


0.132 


0.214 


-0.090 


0.174 


-0.066 


0.143 


15 


-0.037 


0.069 


-0.021 


0.046 


0.124 


0.255 


0.150 


0.231 


-0.094 


0.177 


-0.064 


0.143 



*n = lOr cases are observed by RSS scheme with 10 samples and r replication. 



where 

f 1 + if 1 < i < to 

d.j=< 2 iiTO+l<i</c — to . 

I 1 + if fc - m + 1< i < fc 

Table[T]shows the values of simulated Minimum RMSE (MRMSE) and Minimum Absolute Bias 
(MAB) of H} nn and H^ nn and optimal to for k = 10 and three famous distributions with different 
values of r. From this values one can conclude that Hl nn is better estimator in the sense of RMSE 
and bias. Tables [2] and [3] show the values of Monte Carlo biases and RMSE for H} nn in three 
distributions for n = 10, 20 and 30. This values present a distinct improvement of the estimator 
in RSS scheme relative to SRS scheme. 



3 Goodness of fit tests 

Park, S. and D. (2003) derived the nonparametric distribution function of H c (n,m) as 



„{x) 



if x < rjx or x > rj n + 



-i i 



if m < x < r] i+1 ,i = 1, ...,n ' 



where 



Vi = { ^l( x (i-m) + ••• +»(i+m-l)) if to+ 1 < « < n - to + 1 

+ Efe=n-m+2 n+m -k+l ( X (n) ~ ^(fc-m-1)) X n - m + 2 < i < n + 1 

They used it to correct the moments of the distribution which are used in goodness of fit tests. 

In the exponentiality test, the aforementioned nonparametric distribution is used to estimate 
the mean and A c . 

/(<?:/)= !°° g{x)\n°P-dx (8) 

T c = 1 + log \ C -H c (n,m) (9) 

The following alternatives of the exponentiality null hypothesis have been considered to estimate 
the powers. 

1. Gamma distribution with pdf 

x a ~ x exp(— x) 

f(x;a) = t-V a>0,a;>0 (10) 

7(a) 

2. Weibull distribution with pdf 

f(x; 13) = Px^ 1 exp(-a; /3 ) f3 > 0, x > (11) 



(» 



3. Log-normal distribution with pdl 

1 1 2 

f(x;a) = — -=^exp(- — (logic) ) fj>0,a;>0 (12) 

4. Uniform distribution with pdf 

f(x) = 1 < x < 1 (13) 

As mentioned by Arizono and Ohta (1989), an estimate for /(/, /o), when /o is the normal pdf 
with known parameters /i and a is obtained as 



J mn =log(VW) + -) -H(n,m). (14) 



2n ^—f V o 



2 



When both /j, and a are unknown, we place their estimates, that is, fi = X and a = — X)™=i 
X) 2 in (|29l) and derive the test statistic as 



T = log(V27rd- 2 ) + 0.5 - H(n, m) (15) 

Park, S. and D. replaced the estimates H (n, m) and a with their corrected estimators H c (n, m) 
and a r and derived the test statistic 



T c = log( v / 2^I) + 0.5- H c (n,m) (16) 
In the normality test the following alternatives are considered to estimate the powers 

1. Uniform distribution with pdf 

f(x) = 1 < x < 1 (17) 

2. Chi-square distribution with pdf 

f{x]a) = f(^) { l r/2x{a/2) ' leM ~l x) a>0 > x>0 ( 18 ) 

3. t-student distribution with pdf 

/(a;y)= r( r(y 2 r^u (i+^) ( - +i)/2 » >2 >~ oo<x<o0 ™ 

4. Exponential distribution with pdf 

f{x;X) = Aexp(-Ax), A > 0, x > 0. (20) 



7 



Tabic 4: Critical values for different values of n, m and a 
Exponentiality Normality 



n 


m 








a 










a 




0.1 


0.05 


0.025 


0.01 




0.1 


0.05 


0.025 


0.01 


10 


1 





5357 





6318 





7297 





8617 





5898 





7027 





8034 





9215 




2 


o 


2898 





3546 





4213 





5099 





3765 





4404 





5113 





6005 




3 





2095 





2645 





3243 





3944 





3214 





3712 





4182 





4667 




4 





1619 





2154 





2596 





3293 





3001 





3221 





3593 





3987 




5 





1416 





1916 





2487 





3122 





2903 





3091 





3311 





3544 


20* 


1 





4455 





5091 





5695 





6373 





4775 





5405 





6025 





6587 




2 





2391 





2822 





3305 





3813 





2824 





3264 





3621 





4092 




3 





1707 





2089 





2450 





2939 





2296 





2614 





2940 





3460 




4 





1389 





1738 





2064 





2498 





2073 





2339 





2671 





3112 




5 





1117 





1445 





1772 





2173 





2000 





2287 





2549 





2875 




6 





0964 





1269 





1569 





1918 





1977 





2255 





2501 





2802 




7 





0779 





1114 





1441 





1741 





1968 





2223 





2463 





2695 




8 





0643 





0983 





1250 





1754 





2013 





2225 





2407 





2612 




9 





0495 





0915 





1188 





1604 





2023 





2213 





2375 





2569 




10 





0368 





0797 





1167 





1543 





2008 





2175 





2391 





2512 


30* 


1 





4100 





4567 





4961 





5656 





4273 





4729 





5155 





5768 




2 





2171 





2498 





2796 





3156 





2443 





2776 





3028 





3451 




3 





1516 





1819 





2122 





2402 





1891 





2145 





2408 





2777 




4 





1202 





1481 





1693 





2065 





1649 





1855 





2099 





2407 




5 





0979 





1210 





1503 





1860 





1498 





1739 





1912 





2300 




6 





0825 





1102 





1361 





1559 





1454 





1658 





1879 





2208 




7 





0722 





0950 





1212 





1501 





1433 





1660 





1891 





2134 




8 





0574 





0849 





1061 





1449 





1425 





1663 





1848 





2082 




9 





0574 





0849 





0960 





1270 





1453 





1631 





1833 





2046 




10 





0379 





0635 





0878 





1162 





1428 





1654 





1838 





2039 




11 





0280 





0545 





0760 





1097 





1468 





1661 





1838 





2056 




12 





0151 





0447 





0706 





1030 





1489 





1697 





1875 





2049 




13 





0043 





0354 





0640 





0927 





1502 





1719 





1891 





2072 




14 


-0.0046 





0274 





0612 





0882 





1527 





1720 





1857 





2073 




15 


-0.0190 





0182 





0505 





0813 





1492 





1716 





1871 





2091 



cases are observed by RSS scheme with 10 samples and r replication. 
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Stokes (1980) proposed the sample variance as an estimator of the population variance as follows 



»=i i=i 



This estimator is asymptotically unbiased and asymptotically more efficient than the sample vari- 
ance in SRS. MacEachern et al. (2002) proposed an unbiased estimator of the variance as follows 

a 2 = — (k — 1)MST + (rk — k + 1)MSE, (22) 
rk 



where 

MST 



E E ( x m -Af-^EE i x m - %.) 2 ' (23) 

i 3 3 i 

^EE(%-%.) 2 . ( 24 ) 

3 « 

X[3].=J2 X y]^ r - (25) 



MSE = H r-1) 



and 



A = EE%/ ri (26) 

» j 

If we use our entropy estimator for estimation of Kullback-Leibler distance between an unknown 
pdf and the pdf of the normal distribution, we derive 

7, m „=log(VW) + i-£(^-^J -H 2 mn . (27) 



In goodness of fit test of normality when fj, and er are unknown we can place their estimators in 
the RSS scheme, i.e. /t in (|2"6"|) and the Stokes estimator, (|2T|) to derive the test statistic as 



^ n = log(V2^)+0.5- J ff J 2 nn . (28) 



If we place the MacEachern et al. estimator of variance in ([27]) , we derive another test statistic as 

-, n / _ - \ 2 

KLl n = log(V^) + ±^(^) -H 2 mn . (29) 

Table 2] contains critical values of exponentiality and normality tests for different values of n, 
m and a. 

Table [5] propose a comparison of powers in RSS and SRS schemes, for exponentiality and 
normality tests. The SRS values of powers are given from Park. S. and D. with the modified 
sample entropy of Ebrahimi et al. and their modified estimators of moments. We used the similar 
window size m for the comparison although our maximum powers may be obtained for different 
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Table 5: Power comparison of 0.05 tests against some alternatives in SRS and RSS schemes 

Exponcntiality 



n 







20* 






50* 








(m = 4) 






(m = 6) 




Alternatives 


SRS 




RSS 


SRS 




RSS 


Gamma (1.5) 


0.2176 




0.2740 


0.3480 




0.4193 


Lognormal (1) 


0.2685 




0.1908 


0.6613 




0.4156 


Weibull (1.5) 


0.4639 




0.6199 


0.7752 




0.9056 


Gamma (2) 


0.4862 




0.6218 


0.8281 




0.9050 


Gamma (3) 


0.8816 




0.9693 


0.9993 




0.9999 


Uniform 


0.8021 




0.9979 


0.9989 




1.0000 


Weibull (2) 


0.9138 




0.9896 


0.9995 




1.0000 


Lognormal (0.5) 


0.9967 




0.9994 


1.0000 




1.0000 


Average power 


0.6288 




0.7078 


0.8263 




0.8307 


Normality 








n 












20* 






50* 








(m = 3) 






(m = 4) 




Alternatives 


SRS 


RSS 




SRS 


RSS 










KT 1 




KT 2 


KT 1 


t(5) 


0.1069 


0.0847 


0.0865 


0.2395 


0.1515 


0.1497 


t(3) 


0.1989 


0.1761 


0.1748 


0.5132 


0.4009 


0.3868 


Uniform 


0.3851 


0.4801 


0.4897 


0.8850 


0.9843 


0.9800 


xl 


0.5058 


0.5704 


0.5739 


0.9326 


0.9709 


0.9710 


\\ (Exponential) 


0.8656 


0.9574 


0.9650 


0.9997 


1.0000 


1.0000 


Xl 


0.9934 


0.9999 


0.9999 


1.0000 


1.0000 


1.0000 


Average power 


0.5093 


0.5448 


0.5483 


0.5713 


0.7513 


0.7479 



n = lOr cases are observed by RSS scheme with 10 samples and r replication. 
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Table 6: Maximum powers (maximal m) of 0.05 tests against some alternatives of the null hypothesis 
distributions 

Exponentiality 
n 





10 


20* 


30* 


40* 


50* 


Alternatives 












Gamma (1.5) 


0.5760(5) 


0.3761(8) 


0.4126(12) 


0.4371(14) 


0.4569(7) 


Lognormal (1) 


0.1333(2) 


0.2140(3) 


0.3133(3) 


0.4143(4) 


0.5174(3) 


Weibull (1.5) 


0.5999(5) 


0.7600(8) 


0.8365(15) 


0.8725(8) 


0.9099(7) 


Gamma (2) 


0.5638(4) 


0.7216(8) 


0.7939(7) 


0.8570(8) 


0.9153(7) 


Gamma (3) 


0.9023(5) 


0.9745(5) 


0.9956(5) 


0.9997(5) 


1.0000(3-5) 


Uniform 


0.9201(5) 


1.0000(8,10) 


1.0000(3-15) 


1.0000(2-20) 


1.0000(2-25) 


Weibull (2) 


0.9659(5) 


0.9963(8) 


0.9998(8,12) 


1.0000(4-12) 


1.0000(3-12) 


Lognormal (0.5) 


0.9815(4) 


0.9995(3) 


1.0000(2-6) 


1.0000(2-9) 


1.0000(2-12) 


Normality 








;? 








10 


20* 


30* 


40* 


50* 


Alternatives 












t(5) 


0.0813(4) 


0.0865(3) 


0.1133(3) 


0.1427(2) 


0.1615(2) 


t(3) 


0.1335(4) 


0.1846(2) 


0.2820(2) 


0.3514(3) 


0.4260(3) 


Uniform 


0.1523(2) 


0.5805(10) 


0.9036(11) 


0.9875(16) 


0.9992(16,20) 


xl 


0.3462(4) 


0.6164(4) 


0.8305(6) 


0.9334(6) 


0.9781(5) 


x\ (Exponential) 


0.6926(4) 


0.9670(4) 


0.9992(4) 


1.0000(3-8) 


1.0000(1-14) 


xl 


0.9492(3) 


0.9999(3-5) 


1.0000(1-12) 


1.0000(1-17) 


1.0000(1-22) 



n = lOr cases are observed by RSS scheme with 10 samples and r replication. 



values of m. For normality test two test statistics KL x mn and KL^ nn are compared in the sense 
of power. For n — 20, using the statistic KL? mn cause less powers than KL x mn . Although the 
average of powers of KL? mn gets larger than the average power of KL] nn when n increases to 50, 
but the difference between this powers is ignorable. Since obtaining the statistic KL^ nn is more 
complicated than KLl nn , we prefer to use KL\ nn for the remaining of the study. 

Table |6] shows the maximum powers and the maximal window size, m for a = 0.05 of exponen- 
tiality and normality tests. Ebrahimi et al. (1992) used such maximality to obtain some optimal 
window size m for each n. Table |6] shows that here this values of optimal m differs distinctly 
for different alternatives. In fact choosing an optimal m depends very closely to the alternative 
which is unknown. So in this paper we use the average of powers for considered alternatives as a 
measure to decide about the optimal m. The values of average powers are tabulated in Table [7l 
The authors believe that this values are more useful for the experimenter who wants to perform a 
test, since he is not aware about the alternative. Table [5] shows the optimal m and the maximum 
average powers for different values of n of exponentiality and normality tests. 
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Table 7: Average powers a = 0.05 for different alternatives and different values of n and m 



Exponcntiality 



n 


m 


AP 


n 


m 


AP 


n 


m 


AP 


n 


m 


AP 


n 


m 


AP 


10 


1 


0.2905 


30* 


1 


0.5470 


40* 


1 


0.6156 


40* 


16 


0.7692 


50* 


11 


0.7997 




2 


0.5138 




2 


0.6851 




2 


0.7325 




17 


0.7643 




12 


0.7950 




3 


0.6281 




3 


0.7364 




3 


0.7786 




18 


0.7634 




13 


0.7872 




4 


0.6939 




4 


0.7564 




4 


0.8026 




19 


0.7647 




14 


0.7866 




5 


0.7009 




5 


0.7759 




5 


0.8067 




20 


0.7551 




15 


0.7836 


20* 


1 


0.4477 




6 


0.7640 




6 


0.8056 


50* 


1 


0.6509 




16 


0.7802 




2 


0.6245 




7 


0.7685 




7 


0.7972 




2 


0.7628 




17 


0.7772 




3 


0.6845 




8 


0.7630 




8 


0.7970 




3 


0.8207 




18 


0.7749 




4 


0.7078 




9 


0.7447 




9 


0.7872 




4 


0.8334 




19 


0.7755 




5 


0.7277 




10 


0.7634 




10 


0.7838 




5 


0.8392 




20 


0.7725 




6 


0.7342 




11 


0.7598 




11 


0.7786 




6 


0.8307 




21 


0.7696 




7 


0.7382 




12 


0.7591 




12 


0.7732 




7 


0.8308 




22 


0.7672 




8 


0.7406 




13 


0.7553 




13 


0.7766 




8 


0.8184 




23 


0.7722 




9 


0.7321 




14 


0.7521 




14 


0.7732 




9 


0.8142 




24 


0.7709 




10 


0.7259 




15 


0.7504 




15 


0.7662 




10 


0.8042 




25 


0.7660 


Normality 


n 


m 


AP 


n 


m 


AP 


n 


m 


AP 


n 


m 


AP 


n 


m 


AP 


10 


1 


0.2765 


30* 


1 


0.5229 


40* 


1 


0.5776 


40* 


16 


0.6098 


50* 


11 


0.6746 




2 


0.3470 




2 


0.6154 




2 


0.6811 




17 


0.5960 




12 


0.6702 




3 


0.3622 




3 


0.6547 




3 


0.7130 




18 


0.5808 




13 


0.6647 




4 


0.3876 




4 


0.6628 




4 


0.7173 




19 


0.5702 




14 


0.6583 




5 


0.3520 




5 


0.6559 




5 


0.7072 




20 


0.5540 




15 


0.6517 


20* 


1 


0.4276 




6 


0.6518 




6 


0.6995 


50* 


1 


0.6375 




16 


0.6457 




2 


0.5141 




7 


0.6361 




7 


0.6911 




2 


0.7269 




17 


0.6385 




3 


0.5483 




8 


0.6219 




8 


0.6786 




3 


0.7482 




18 


0.6325 




4 


0.5586 




9 


0.6226 




9 


0.6689 




4 


0.7479 




19 


0.6227 




5 


0.5418 




10 


0.6104 




10 


0.6598 




5 


0.7415 




20 


0.6160 




6 


0.5294 




11 


0.5983 




11 


0.6521 




6 


0.7283 




21 


0.6090 




7 


0.5230 




12 


0.5799 




12 


0.6425 




7 


0.7148 




22 


0.5952 




8 


0.5078 




13 


0.5631 




13 


0.6363 




8 


0.7043 




23 


0.5858 




9 


0.4891 




14 


0.5505 




14 


0.6258 




9 


0.6903 




24 


0.5713 




10 


0.4744 




15 


0.5291 




15 


0.6168 




10 


0.6827 




25 


0.5592 



n = lOr cases are observed by RSS scheme with 10 samples and r replication. 
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Table 8: Values of the window size m with largest average of powers against alternatives 



n 


Optimal m(max average power) 




Exponcntiality 


Normality 


10 


5(0.7009) 


4(0.3876) 


20 


8(0.7406) 


4(0.5586) 


30 


5(0.7759) 


4(0.6628) 


40 


5(0.8067) 


4(0.7173) 


50 


5(0.8392) 


3(0.7482) 
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